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Abstract 

We transform the quartic Hubbard terms in the extended Hubbard model to a 
quadratic form by making the Hubbard-Stratonovich transformation for the electron 
operators. This transformation allows us to derive exact results for mass operator 
and charge-charge and spin-spin correlation functions for s-wave superconductivity. 
We discuss the application of the method to the d-wave superconductivity. 
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1 Introduction 



The Hubbard model predicts pliase instabilities which give rise to a divergence 
of the charge and spin correlation functions, and therefore, it has been the 
focus of particular interest as a model for high-temperature superconductivity. 
The Hamiltonian of the standard Hubbard model contains only two terms 
representing the hopping of electrons between sites of the lattice and their 
on-site interaction. If the interaction between electrons on different sites of 
the lattice is included, the model is referred as the extended Hubbard model. 
In what follows we study the following Hamiltonian: 

H = - ^iji'lcri^j,^ - fJ'Yl ^*.<^ + ^ H ^ht^hi ~^ H ni^afijy, (1) 
i,j,<T i,cr i <i,j>a<T' 



Email address: Zlatko.Koinov@utsa.edu (Z.G. Koinov). 



Preprint submitted to Elsevier 



1 February 2010 



where /i is the chemical potential. The Fermi operator -ipj^ i'4'i,a-) creates (de- 
stroys) a fermion on the lattice site i with spin projection a =f , J, along a 
specified direction, and nj^o- = '^lai'i^o- is the density operator on site i. The 
symbol I]<ij> means sum over nearest-neighbor sites. The first term in ([T]) 
is the usual kinetic energy term in a tight-binding approximation, where tij 
is the single electron hopping integral. Depending on the sign of f/, the third 
term describes the on-site repulsive or attractive interaction between electrons 
with opposite spins. We assume that > 0, so the last term is expect to sta- 
bilize the pairing by bringing in a nearest-neighbor attractive interaction. The 
lattice spacing is assumed to be a = 1 and the total number of sites is A^. 
The simplest method to study the possibility for the extended Hubbard model 
to show a superconducting instability is to apply mean-field analysis of pairing 
followed by general random phase approximation (GRPA) [TP] . Going beyond 
the GRPA requires reliable approximation schemes to handle self-consistent 
relations between single- and two-particle quantities: the mass operator E 
depends on the two-particle Green function K, and the kernel of the Bethe- 
Salpeter (BS) equation 5T,/5G for the spectrum of the collective excitations 
itself does depend on the mass operator. A possible approximation to this 
problem is the so-called two-particle self-consistent (TPSC) approach [3]-P^. 
The TPSC approach is a method for closing the set of equations for single- 
particle mass operator and the two-body density matrix operator. The later 
can be factorized by introducing the so-called equal-time pair-correlation func- 
tion gaa'ihj) [I^|T3] which itself depends on the density-density correlation 
function. In other words, the TPSC approach goes beyond the GRPA for 
single-particle mass operator by establishing a self-consistency relation be- 
tween single-particle and two-particle quantities. By setting gaa'ihj) = 1 one 
should recover the GRPA results for the mass operator and charge and spin 
correlation functions. 

In what follows, we first obtain exact formulas for the electron self-energy 
(electron mass operator), the charge and spin correlation functions. We also 
briefiy discuss how our approach could be generalized in order to include d- 
wave instabilities of the types examined in Refs. p^l7|18|19|20] . 



2 Field-theoretical approach to extended Hubbard model 



The interaction part of the Hamiltonian ([T]) is quartic in the Grassmann 
fermion fields so the functional integrals cannot be evaluated exactly. However, 
it is convenient to transform the quartic Hubbard terms in ([1]) to a quadratic 
form by making the Hubbard-Stratonovich transformation for the electron 
operators: 
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fi[A] exp 



The symbol "hat" over any quantity O means that this quantity is a matrix. 
The functional measure -D/u[A] is chosen to be: 

The Hubbard- Stratonovich transformation converts the quartic problem of 
interacting electrons to the more tractable quadratic problem of noninteracting 
Nambu fermion fields 




coupled to a Bose field Aa{z) where a =t, i is the spin degree of freedom which 
reflects the spin-dependent nature of the Hubbard interaction. The bare boson 
propagator in ([21) provides an instantaneous spin-dependent interaction, and 
in accordance with the Hamiltonian ([1]), it should have the following form: 



D^^liz, z') = D%,j'; v-v') = S{v - v') [US,,,S^-^ - 2l^<,,,> (5^,^ + <5„^) 



1 ^^e{4k-(r."r,)-..(W)]}^(o)(k.^^^)^ 



Di°j(k; tUp) = US^^^ - V{k) + 5a,/3) . 

Here a is complimentary of a, and V^(k) = AV (cos + cos ky) is the nearest- 
neighbor interaction in momentum space. The symbol V^<jj'> is equal to V if j 
and j' sites are nearest neighbors, and zero otherwise. We have used composite 
variables y = {ri,u} = {i,u}, x = {ri',u'} = {i',u'}, z = {rj,v} = {j,v} and 
z' = {rj',v'} = {j',v'}, where rj,rj/,rj and rj/ are the lattice site vectors. 
The symbol J2lo is used to denote /3~^ J2p- For boson fields we have Up = 
(27r//3)p;p = 0,±l,±2,.... 

After performing the Hubbard-Stratonovich transformation, the action of the 
system becomes 

5 = 4^) + 4^) + 5(^-^), (4) 



where: 

St^ = ti{y)G^'^-\y-x)i^{x), (5) 
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(6) 
(7) 



The inverse Green function of free electrons G*^"-* ^{y]x) is diagonal with re- 
spect to the spin indices and has its usual form: 



GW-Ht,i/;t,^) 



(8) 



where Gj°j""^(k, lojm) = [loOm - (e(k) - /x)]"^ , and G|°j""^(k, tUm) = [ioom + (e(k) - 
Here e(k) = —2t{coskx + cos ky) — At' cos cos ky is the non-interacting dis- 
persion on a square lattice, fi is the electron chemical potential, and the 
symbol X^^^ is used to denote (5~^ J2m- For fermion fields we have — 
{27i/P){m+l/2)im = 0,1,2, . 

The bare vertex T^aKVi a; | -z) is a 2 x 2 matrix defined as follows: 







-r(^\i,y;i,x\z)) 

rW((7, y;a,x\z)^ r^") (a, i, u; a, i', u' \ i", v) 
^ 5{u - v)5{u - u')5a,aSi^i'Si^i". 



(9) 



Since the electrons polarize the boson field, and the boson field acts onto the 
electrons, our approach describes the correlated motion of the electrons and 
the surrounding polarization field. 

In field theory the expectation value of a general operator 0{u) is expressed 
as a functional integral over the boson field A and the Grassmann fermion 
fields -0 and -0 



< T^{0{u)) > 
exp 



Z[J, M] 

Ja{z)A^(z) - ^{y)M{y;x)i^{x) 



Dfi[ip,'ip,A]d{u) X 



\J=M=0: 



(10) 
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where the symbol < ... > means that the thermodynamic average is made, and 
Tu is an m— ordering operator. J, M are the sources of the boson and fermion 
fields, respectively. The functional Z[J,M] is defined by 

Z[J, M]= J DIJ§, ^i^^[Mz)A^{z)-^{y)M(y-,x)ii(x)]^ ^^^-^ 



where the functional measure D^l^p^'ip, A\ = DADipBtpexplS) satisfies the 
condition / D fi[tp , tjj , A] = 1. 

It is convenient to introduce complex indices 1 = {(Ti,xi}, 2 = {o"2, 2/2}, 
where, ai^2 = {tii} and Xi = {rjj,Mi}, and j/2 = {ri^,U2}. We define a 
functional derivative 5/5M(l; 2), and depending on the spin degrees of freedom 
0"! and o"2, there are four possible derivatives: 



(5M(;,|/2;t,a;i)' 6M{i,y2;l,xi)' 
The reason to write the expectation value (fTOj) as a functional integral is that 
all Green functions related to system under consideration can be expressed 
in terms of the functional derivatives of the generating functional of the con- 
nected Green functions iy[J, M] = lnZ[J, M]. By means of the functional 
iy[J, M], we define the following Green and vertex functions of the extended 
Hubbard model: 
Boson Green function: 

The single-electron Green function G{1]2) = —6W/6M{2] 1) in the Hubbard 
model assumes the form: 



_ I < r„ f^^(xi)^|(?/2)j > <T^{^t{xi)ilj^{y2)) > 

< r„ (^|(xi)^|(l/2)) >< T„ (^|(xi)V^^(y2)) > 



Depending on the two spin degrees of freedom ai and (T2, there exist two 
"normal" Green functions 

G(t,xi; 1,1/2) = Gt,t(a;i,y2) = - < Njt{xi)^\{y2)) >= -jvttt r r, 



G{l,xi;l,y2) = Gi,i{xi,y2) = - < fu ['ipl{xi)^i{y2)) 



>-- 



6W 



6M{i,y2;i,xi)'' 
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and two "anomalous" Green functions 



6W 

>= 



6W 



1,1/2) = G4,t(a;i,Z/2) = - < f „ (V^J(xi)^j(y2 
G{t,xi;l,y2) = Gi^,l{xi,y2) = -<fuiij^{xi)iJi{y2))> 

OlVl 1/2, I , Xi) 

We introduce Fourier transforms of the "normal" G-^^^lhjU — «') = — < 

^^«(^t,k("')^|,k'^^)) '^u(k,M-^^') = - < f„(^jj^(M')^;,k(")) > 
"anomalous" Gi^-f(k,u — u') = — < Tu{i'^i.{u')ip^]^{u)) >, G^^i(k,u — u') = 
— < k('"')V^j kl"")) > one-particle Green functions. Here ip^-^^iu) , ip^ ]^(n) 

and V'^j^l'"), kl"") ^ire the creation- annihilation Heisenberg operators. The 
final form of single-particle Green function is given by 



G(l;2) 



^^exp{2 [k. (rj^ - r,J -M2)]} 



(14) 



^^^^(k,^^^) G^^i(k,iujm) 
^^^^(k,^^^) G^^i{k,tuj^) 

The two-particle Green function is defined by: 
/ 



K 



1 3 



6'W 



_ _ 6G{1;2) _ 
2 4 I ~ 5M(2;1)5M(3;4) ~ 5M(3;4)' 



(15) 



Depending on the four spin degrees of freedom ai,cr2,o"3 and (J4, there are 
sixteen different components of the two-particle Green function. 
The vertex function ro-^(2; 1 | 2;) is given by: 



r.(2;i|.)^-^gfilz.,-.V,.) 

= G-\2;3)^-^^G-\A;l)D,^{z\z). (16) 

If the spin variable a is fixed, then depending on the spin degrees of freedom 
0"! and (72, there are four different vertex functions. 

The action (j4j) as well as all of the above definitions allow us to map the 
extended Hubbard model onto the polariton model that describes the light 
propagation in semiconductors [21]. This mapping allows us to apply directly 
to the Hubbard model all exact equations and relationships derived for the 
case of light propagation in crystals. For example, we can write the mass 
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operator S(l;2) = S//(l;2) + Ilp{l;2) as a sum of Hartree and Fock parts 
|21j . The Hartree contribution to the mass operator is diagonal with respect 
to the spin indices: 



= rW(l;2 I z)G(4;3)r(°)(3;4 I z')dI^I{z, z') = 5(1-2) x 

{-UG^^^^^{ii,Ui]ii,Ui) 
-2VY^[G^^^^^{ii,Ui]ii + a,ui) - G^^^-^^{ii,Ui]ii + a,Ui)]}, (17) 

a 

where the summation on a runs over the nearest-neighbor sites of site ii. 
The Fock part depends on the boson Green function D, or equivalently, de- 
pends on the two-particle Green function K: 



Sir(l; 2) = Si7'(o-i, 2i, Mi; (72, i2, U2) 

= -ri°)(l;3 I z)G'(3;4)r^(4;2 I z')D^p{z, z'), 

(53' 

-r^^\l;Q\z)D^^I{z,z')rf{A;5\z')K 1 G-^(3;2) 
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(CTi, Zi, Ml (T3, ^3, M3 \ ^ ^ 

' ' 3, 3, 3 I ^2VJ2[K 



(Ti,ii a, Ml 0-3,23, Ms 
cri,ii + a, Ml (71,21, Ml 



^ ai,ii + a. Ml Cr3,i3,M3 



]}G (0-3, 23, M3; 0-2,22, M2)- 



y 0-1, 2i a. Ml 0-1, 2i, Ml 

The \^-terms in (fTSi) differ from the corresponding result in [T2][T3] . 



3 Spectrum of the collective modes 



Spectrum of the two-particle excitations (or collective modes) c<;(Q) can be ob- 
tained by locating the positions of the common poles of the Fourier transform 
of the two-particle fermion Green function fllSp and the Fourier transform of 
the boson Green function fll2p . In other words, collective modes are defined 
by the solutions of the BS equation for the function K or the Dyson equation 
for the boson function D. 



The BS equation is ^\E' 



7^{o)-i _ J 



\E' = 0, where / is the kernel, and 



-1 



1 3 

2 4 



^-^(1; 3)^-^(4; 2 
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is the two-particle free propagator constructed from a pair of fully dressed 
single-particle Green functions. The kernel / = + depends on the 
functional derivative of the Fock contribution to the mass operator. Since 
the Fock term itself depends on the two-particle Green function K (see Eq. 
f|T8l) ). we have to solve self-consistently a set of two equations, namely the 
BS equation and the Dyson equation G^^ = G'^^-'^^ — S for the single-particle 
Green function. 

Similar obstruction arises if we start from the Dyson equation for boson Green 
function: 

D-'^{z,z') = ^^-^(z,/) -n„;3(^,^'), (19) 



where Ilap{z, z') is the proper self-energy of the boson field. The obstruction 
now is that the proper self-energy depends on the vertex function F, or on the 
two-particle Green function K: 



Ii^p{z,z') = V^^\l;2\z)K 



^2 3 



Tf{'i-A\z') 

1 4 " 



= Ff (1;2 I ^)G(2;3)G(4;1)F^(3;4 I z'). (20) 

The Green function K, which we shall call the Green function of electronic 
excitations, satisfies the following BS equation: 



Let us introduce the so-called general response function lias'- 



Ii^p{z-z') = V^^\l-2\z)K\ \vf{?>-A\z'). (22) 





By means of (1221) . we can rewrite the BS and Dyson equations as follows: K = 
K(o) + K(°)F(,°)n„;3F^°^i\:(°), D = + DWnD(°). The Fourier transforms of 
the general response function Ua/siQ; ou) and the proper self-energy I[a/3{Q;uj) 
are connected by the following equation: 



n„/3(Q; uj) = n„;3(Q; oj) + n,^(Q; to) i-U + V^(Q)) %^(Q; u) 

+U^^iQ;uj)V{Q)U^p{Q;u), (23) 
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According to Eq. ( 120|) . the proper self-energy 11, the Green function K and 
the vertex function F must have common poles. Let -^/(Q) and Q denote the 
energy and momentum of these common poles. Close to -E'«(Q) one can write: 













K 




= K 




^2a) 






i2,U2 



-ZUJp{ui-U3) 



(24) 



where $o-T;o-2(ri2, Ui — U2) are the BS amplitudes: 



iQ 



[ri^,ri^;ui - U2) = exp 



Due to the form of the bare vertex f "^^^ we have to take into account only the 
equal "time" ui = U2 amplitudes: 

0ji,^(r,, - r,2; 0) = 1 J^expjzk. (r,, - v,,)}<Pi^^^^{]^, Q), 



where 0j^-f-(k, Q), ^^-^^(k, Q), 0j^(k, Q) and (t)\i{}^, Q) are the equal "time" two- 
particle wave functions in k-representation. By means of ( 12^ and ( l20l) we 
obtain: 



na/3(Q,Cu) = 



0S(O; O)0;3*(O; 0) 0S*(O; 0)0^(0; 0) 



(25) 



where w^q = Ei{Q) — /i. If the proper self-energy is known, the spectrum of 
the collective excitations w(Q) could be obtained assuming the vanishing of 
the following 2x2 determinant: 

det II 5^,(s - i-U + ViQ)) U^^iQ, u) + V(Q)n„^(Q, u) \\= 0. (26) 



By solving Eq. fl26l) (with the help of Ua/s = H^q) we find two different types 
of collective modes. The first one is governed by U interaction: 



= 1 - f/ 



nn(Q,a;)-nt4(Q,a;) . (27) 



Strictly speaking the V interaction is included indirectly in Eq. (127|) through 
BS amplitudes and poles of the function K. This collective mode manifests it- 
self as a pole of the spin response function Xss{Q] u) = 2 [Il^^{Q] uj) — n-f-4,(Q; u)], 
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and therefore, gives rise to the spin instabihties. The spin correlation function 
is defined by 

Xss{,Q]Uj) = - ^ , . . , 28 

1 - Xs4Q; w)f//2 



where Xss{Q; w) = 2 ntt(Q; - n^lQ; 

The second type of collective modes satisfies the following equation: 



= 1 + [U- 2Vm U^iQ, cu) + U^^iQ, CO] 



(29) 



The second collective mode manifests itself as a pole of the charge response 
function Xcc{Q',<^) = 2 [n^(Q; w) + n|4.(Q; w)] which is determined by the 
following equations: 



Xcc{Q;uj) 



i + XcciQ;uj) [U/2-vm' 



(30) 



where Xcc(Q; w) = 2 [n^(Q; w) + UniQ'^ ^) 
Within the GRPA, one should replace K in ( !20l) by K^^\ thus obtaining the 
free response functions instead of the exact x expressions. In other words, 
the exact relations ( l28l) and (130|) in the GRPA are given by 



Xss{Q;uj) 
XcciQ;^) 



xf;(Q;^) 



xi?(Q;a;)f//2' 

x(?(Q;c^) 



l + xic^(Q;u;) [U/2-Vm 



(31) 
(32) 



4 Discussion 



We have established exact results for mass operator and spin and charge cor- 
relation functions in the case of s-wave superconductivity described by the 
extended Hubbard Hamiltonian. The Hamiltonian depends on an on-site re- 
pulsive interactions U, which drives antiferromagnetism, and a near-neighbor 
attractive interactions V, which drives d-wave superconductivity. Let us dis- 
cuss consider the extended Hubbard model with a nearest-neighbor repulsion 
V. In this case the near-neighbor repulsive interaction V drives instabilities 
related to the change in Fermi surface topology (Pomeranchuk instability) 
[TSflQ] . The Fermi surface is related to the diagonal elements of matrix (|T8l) . 
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In the RPA we replace K \ by ^(1; 3)^(4; 2), and therefore, the U 

term in (JTSl) does not create any changes in the non-interacting Fermi sur- 
face /i = e(k). In RPA only the near-neighbor repulsive interaction (V— ?■ -V) 
changes the Fermi surface topology through the diagonal elements of the mass 
operator (|T8l) : 

(Ej7(k))^^ = -AVj2[cosih - Qx) + cos{ky - qy)] n^(q), (33) 

q 



where na(k) = Y.ujm ^^^(k, ioom) is the occupation of the site k in momentum 
space. Note that the V^-terms of Ej?(l,2) in Refs. [T^P3] are proportional to 
the delta function 5(1 — 2), and therefore, Ei?(k) does not depend on k. 

It is found [TS] that the near-neighbor repulsive interaction V enhance small 
anisotropies producing deformations of the Fermi surface that break the point 
group symmetry of the square lattice at the Van Hove filling. Since the antifer- 
romagnetic order is favored by U but suppressed by V, one has to expect that 
the Pomeranchuk instability competes with magnetic instabilities and will be 
suppressed at some critical value V. 

The analog of the BCS reduced Hamiltonian for singlet density-wave order is 
the so-called f-Hamiltonian [T71I19II2UI: 



^ = II(e(k) - + — S{q)di.dpn^{q)np{-q), (34) 

k,cr k,p,q 



where (ij^ = cos{kx) — cos{ky) and the operator n]^(q) is n]^(q) = J^ai^ti^ + 
|q)'?/'o-(k — |q). In momentum space the corresponding order parameter can 
be defined through the following operator (for more general definition of 0(q), 
see Ref. [20]) 

0(q) = ^dkV^t(k+ lq)5^^(k_ Iq)^ (35) 
k 



where ^/'(k) and V'(k) are the Fourier transforms of the Numbu fermion fields, 
and is the Pauli matrix. The interaction term in the above Hamiltonian in 
coordinate space is ^J2i.jS{ri — rj)0{ri)0{Tj), where the Fourier transform 
of S{ri) is S'(k). Now, we introduce a boson field B{z) = B{rj;v) with a 
free propagator S^^\z, z') = S{v — v')S^^\rj — vji) and a bare vertex function 
T^'^\y,x\z) = T^^\ri,u;ri^,u'\rj,v) = 5{u — u')S{u — v)S{ri — rii)5{ri — rj). Now, 
we can perform exactly the same steps as in the case of s-wave superconduc- 
tivity. As a result, one can obtain the following equation for d-wave response 
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function n(Q, u): 

n(Q, u) = n(Q, u) + n(Q, a;)5(°nQ)n(Q, w), 



(36) 



where the proper self-energy is defined as 

n(z, z') = r(°)(l; 2\z)G{2; 3)G(3; 4)1(3; 4\z). 
Solving ( |36|) we obtain the GRPA result for the d-wave response function: 

(n ^ x^'KQ,^) 

X[^^^) i_5(0)(Q);^(0)(Q,,,)' 
where is the d-wave proper self-energy in RPA (see eq. (7) in Ref. [19]). 

In summary, we have used the Hubbard-Stratonovich transformation to con- 
vert the quartic Hubbard problem of interacting electrons into more tractable 
quadratic problem of noninteracting electrons coupled to a Bose field. This 
field-theoretical approach allows us to express all quantities of interest in 
terms of the corresponding Green functions and obtain exact relations be- 
tween single- and two-particle quantities. 
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